Growth of aqueous foam on flexible membranes 
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In this paper, I study the coarsening dynamics of two-dimensional dry foam sandwiched by de- 
formable membranes. The time- varying deformation of the confining membranes gives rise to a 
significant alteration in the evolution of polygonal cells of bubbles when compared to the case of 
rigid membranes. This alteration is attributed to the correlation between the rate of inter-cell gas 
transfer and temporal fluctuation in surface curvature within a cell domain. The existing material 
constants are referred to understand the utility of the correlation effect toward the artificial control 
of the coarsening dynamics. 
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I. INTRODUCTION 

Aqueous foam exhibits a good interplay between geom- 
etry and physics. With time, foam consisting of polyhe- 
dral bubbles evolves into the equiUbrium structure, dur- 
ing which internal gas diffuses from a bubble to others 
through thin curved liquid interfaces [H, [3j S 3 • Diffu- 
sion is driven by the pressure difference between bubbles; 
assuming the constant diffusion coefficient, the pressure 
difference in two adjacent bubbles is proportional to the 
geometric curvature of their common boundary interface. 
Each boundary moves toward its concave sides due to 
the inter-bubble gas transfer, where the velocity of the 
boundary motion is again proportional to curvature Q. 
As a result, some bubbles dilate while others shrink and 
eventually disappear, which results in a progressive in- 
crease in the average bubble size, i.e., the coarsening of 
foam. 

Foam that we encounter in our daily life, such as shav- 
ing cream and beer head, consists of a three-dimensional 
agglomerate. Its coarsening dynamics [1, 0, H H [13, [lHj 
[l2, 13.,Ii,i5.,Ifi,iIJ8Ual2p,[2^ as well as equihbrium 



cellular structures [22, 123, |2J| have required formidable 
efforts for clarification because of geometrical and topo- 
logical complexity. This is partly the reason why a large 
degree of attention has been paid to two-dimensional 
counterparts, i.e., a foam monolayer confined between 
two membranes ^ [H, [13, H, ^ ^ [HI, [13, [H, [H. 
It is remarkable that in two-dimensional foam confined 
between flat planes, the time evolution of a polygonal 
cell of bubble depends only on the number of its sides, 
regardless of its shape or area. The growth-rate of the 
area S of an ri-sided cell is given by dS/dt oc (n — 6) 
[35| . This formula states that a cell is stationary if n = 6 
but it grows (shrinks) at a constant rate if n is larger 
(smaller) than 6. The disappearance of shrinking cells 
causes a topological change in the network of liquid in- 
terfaces, whose effects on the stability of foam has also 
been largely investigated [13, [H, [13, [11] . 
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The coarsening behavior on the flat plane alters dras- 
tically when the foam is constrained to a curved surface 
[sol. [40|. In the latter case, the evolution of cells is char- 
acterized by the Gaussian curvature K of the underlying 
surface. When the surface has a positive (negative) cur- 
vature, 71-sided cells with n < 6 (ti > 6) can be stationary, 
yielding dS/dt = 0, only if S equals to a specific value 
that depends on n and K (see Eq. ((5|) below). Further- 
more, the stability of those stationary cells is sensitive to 
the sign of K . For instance, no cell on a positively curved 
surface is stable; once a cell grows (shrinks) slightly un- 
der perturbation, then it keeps growing (shrinking). In 
contrast, all stationary cells on a negatively curved sur- 
face are stable; therefore, the equilibrium configuration 
consists of various n-sided cells each having a specific 
area determined by n(> 6) and K. Such two-dimensional 
foam spreading over a curved surface could be realized on 
an elastic confining plate or on a phase boundary with 
another fluid medium that repels the foam. 

The present article provides a further generalization of 
the coarsening on curved surfaces, i.e., the coarsening of 
foam constrained on a flexible surface exhibiting time- 
varying deformation. The successive deformation of the 
confining surface changes the pressure difference of adja- 
cent cells that drives gas transfer across liquid interfaces. 
As a consequence of the correlation between surface de- 
formation and inter-cell gas transfer, the growth rate of 
the cell shows an intrinsic difference from that of the 
rigid curved surface. Realistic material parameters are 
employed to prove that the correlation effect plays a dom- 
inant role in the actual coarsening dynamics on flexible 
membranes under appropriate physical conditions. 



II. COARSENING ON A RIGID CURVED 
MEMBRANE 

This section briefly reviews the coarsening theory of 
two-dimensional dry foam on a rigid curved surface. The 
term "dry" refers to the assumption that liquid films be- 
tween adjacent cells are so thin that they can be treated 
as curves with no thickness and the vertices can be 
treated as points. This assumption has long succeeded 
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in exploring the nature of coarsening dynamics of foam, 
while considering that the realistic shapes of films with 
finite thickness [4ll, 13] and their effects on permi- 

ability [i^ may encourage quantitative agreements with 
experiments; we shall revisit this point in Section V. 

Let us assume that the monolayer foam is confined be- 
tween two rigid membranes with spatially uniform Gaus- 
sian curvature K. The gap between the membranes is 
smaller than the typical length of boundary curves. The 
growth-rate of the area of an n-sided cell with internal 
pressure p is described by 



d5 

dt 



(1) 



where Apj = p — pj is the pressure difference between 
the cell and its jth neighbor, £j is the length of the jth 
boundary curve separating the two cells, and 7 > is 
a diffusion constant. Equation (jT]) captures the simple 
idea that if the cell has a higher pressure than the jth 
neighbor {i.e., p > Pj), then gas escapes to the neighbor, 
and vice versa. With local equilibrium, Apj is balanced 
by the line tension a along the interfaces, satisfying the 
generalized Laplace- Young law [3§| 



Apj 



(2) 



where kj is the geodesic curvature of the jth interface. 
From viewpoints of differential geometry, any n-sided 
polygon on a surface with curvature K satisfies Gauss- 
Bonnet's theorem [4^[43| expressed by 



i=l 



{tt- 



3 = i 



KdS = 27r, 



(3) 



where Ui is the internal angle of the ith vertex and must 
be equal to 27r/3 for all i according to Plateau's lemma 
[4^. Consequently, we obtain the result 



dS 



-{n-6) + KS 



(4) 



which is called a generalized von-Ncumann formula de- 
scribing the foam coarsening on a rigid membrane (soj . 

Formula (|4|) accounts for the stability properties of n- 
sided cells on rigid curved surfaces. For K > 0, such cells 
that satisfy n < 6 and 



S = S* 



K 



(5) 



can be stationary, although all stationary cells are un- 
stable. For instance, if S becomes slightly larger than 
S* due to perturbation, then the quantity in the square 
brackets in Eq. (|4|) becomes positive. Therefore, we ob- 
tain dS/dt > after the perturbation, which signifies a 
persistent growth in the cell. On the contrary, all sta- 
tionary cells are stable for K < 0, since S being larger 
(smaller) than S* makes dS/dt negative (positive). In 
this context, the case of a flat plane is marginal, in which 
the stationary cell of n = 6 is neither stable nor unstable 
against perturbation. 



III. COARSENING ON A DEFORMABLE 
MEMBRANE 

Now, we focus our attention on the case where the con- 
fining membranes are mechanically flexible. Membrane 
deformation induces a change in the value of K within 
a cell domain, thereby changing its area by incrementing 
dScur- In addition, inter-cell gas transfer contributes to 
the change in the area by incrementing SSdn, which is 
similar to the case of rigid curved membranes discussed 
in Section II. An important consequence of membrane 
deformation is that it alters the internal pressure of the 
cell, which causes a change in the pressure difference be- 
tween cells that works as a driving force to yield ^S'dif. 
As a result, the membrane deformation and inter-cell gas 
transfer correlate with each other, which imply a sizeable 
deviation in the growth rate equation from Eq. (|4]) that 
we previously derived. The quantitative determination 
of such a correlation effect on the cell growth equation is 
the main purpose of the present work. 

Let us assume that the Gaussian curvature K of the 
confining membrane is spatially uniform and varies con- 
tinuously from K = Kq at t = to i^i = Kq + SK at 
t = St. The increment SS of the cell area obtained at 
t — St is the sum of the two contributions: 



SS = SS, 



dif 



SS, 



where 



and 



SS, 



dSdii 



dif 



dt 



St, 



SS, 



dSc 



dK 



5K, 



(6) 



(7) 



(8) 



with SK = [dK / dt)St. The subscripts and Kq im- 
ply differentiations at i = and K = Kq, respectively. 
dSd\f/dt\Q in Eq. ^ is identified with dS/dt appearing 
in Eq. ([4]), since both describe the growth rate before the 
deformation occurs. Equation ([8]) has an explicit form 
that is given as 



SSc,-\\Y — 



02 

Op 



2 (27r - KqSg 



-SK, 



(9) 



and it shall be proved in Appendix A. From Eqs. ([6|)-([9|), 
the area Si ^ Sq + SS aX t = St is readily evaluated. 

Next, we consider the area growth within the time du- 
ration [(5i, 25^]. The increment SScmv in this interval is 
obtained by replacing Kq and in Eq. © by A'l and 
Si, respectively. On the other hand, some caution is re- 
quired in deriving the form of SSd\i at this stage because 
of the correlation between internal pressure and mem- 
brane deformation. The pressure p of a given cell after 
deformation is expressed by 



Pi = Po + Spdii + Spc 



(10) 
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where Spdn is the pressure increment that would be ob- 
tained provided SK = 0, and Spcur is the one provided 
no diffusion occurs during deformation. Wc denote by 
and the counterparts of the jth neighbor, both of 
which obey the similar expressions of p and S, i.e., 



Pi 



Po 



5p^ 



'dif 



Sp 



cur 7 



S{ 



Sl+5Si^ + 5Si^,. (11) 



To analyze the deformation effect on the pressure differ- 
ence Apj = pi — p\, we take notice of the fact that the 
number of gas molecules in the cell is preserved during 
deformation if no diffusion occurs. This conservation law 
is formally represented by (po+<5pcur)('S'o+<^'S'cm) = PO'5'oj 
or equivalcntly 



Po ■ SScur + SpcuT • 5*0 ~ 0, 



(12) 



which correlates Spcur to SScui. We also sec from the 
Laplace- Young law that 

Po-Po = cri^o, 

{po + Spdii) - {pi + Sp^^if) = <JKl^, (13) 

where Kq is the geodesic curvatures of the jth boundary 
observed at t = 0, and is the fictitious one that would 
be observed ifSK = during [0, St]. From Eqs. pUl) - p^ . 
it follows that 



Ap{ 



"-St 



■PO - — TT- [Po 



So 



(14) 



Substituting Eq. (fTil) in Eq. ([T|) and applying Gauss- 
Bonnet's theorem, wc obtain 



dSdiS 



dt 



St 



7^ [|(n - 6) + KoSo + SKS„ + KoSSdn 



'0 



qj 



Po 



J ^J' 



(15) 



where the second-order terms with respect to increments 
were neglected. It is to be noted that all the increments in 
the right side are those obtained in the previous interval 
[0,5t]. 

Equation (jlSp determines the diffusion-induced incre- 
ment SSdif = dSdif/dt\st5t obtained at t — 26t. The 
area growth for larger t can be evaluated by succes- 
sively applying the procedure shown above. After de- 
ducing 5*2 at t = 2St, for instance, we rewrite the set 
{S2,K2,Si,Ki,pi} by {Si,Ki,So,Ko,po} to calculate 
S2 again, which provides the subsequent value of the area 
at i = 35t. An explicit algorithm of pursuing the time- 
varying p is given in Appendix B. 

It should be emphasized that in Eq. (jTS]) , the deforma- 
tion effect manifests in the product SKSq in the square 
brackets and the two summations with respect to j. In 
particular, the presence of the last two summations in- 
dicates that the diffusion-induced growth rate of a cell 
becomes dependent on the local environment around the 
cell. This situation is in contrast to the case of a rigid 
membrane described by Eq. where the growth rate 
is determined only by the properties of the cell itself. 



IV. DEFORMATION EFFECT ESTIMATION 



To estimate the deformation effect in Eq. (|15p . we re- 
place the fractions 5SI^j./Sq by its mean value over n 
adjacent cells and the sum of edge lengths J^j by the 
perimeter £ of an effective circular domain whose area 
equals to the original polygonal cell area 5*. It follows 
that C Sit t = St is represented by 



(16) 



which will be derived in Appendix A (sec Eq. (|A4p ). 
Straightforward calculation yields 



dSdif 



dt 



7(T 



St 



:{n^6) + KoSo (1-/3) 



-7fT [SKSo + KoSSdif) + jpoi/3 - a)C, (17) 



where 



a 



n 

-E 



3 = 1 



"^0 



(18) 



To proceed with the arguments, we consider sub- 
millimetcr-scale bubbles of S ^ 0.1 mm^ under slightly 
time-varying curvature of SK ^ 10~^ mm~^ per sec- 
ond with K — 1 mm^^ at t = 0; these conditions are 
in the realm of laboratory experiments. Then, we have 
/3 ^ SSi^^./So <C 1, as a result of which Eq. pT]) is sim- 
plified as 



dS, 



dif 



dt 



St 



7(7 

' L3 
lPo{f3 



(n — 6) 
~a)C. 



Kit) (So + SSdii) 



(19) 



The most important deviation of Eq. (fTO|) from Eq. ([3]), 
i.e., from the growth rate equation for rigid membrane 
cases, is the presence of the term jpo{P — a)C. In fact, 
this term relates ^S'cur obtained in the previous time in- 
terval, say, [0, (St], to SSdu obtained in the subsequent 
interval, [St,2St]. When the confining membranes are 
rigid, then this term vanishes since a oc SS cur SK = 0, 
and so does /3. 

The salient finding of Eq. is the fact that by re- 
ferring realistic material constants such as Pq_= 10^ Pa, 
(7 = 10^ N/m, and 7 = 10"^ m/(Pa-sec) 0, Sll, we 
obtain 



and 



IPoiP — ct)C- ^ 10 ^ mm^/sec, 



7cr = 10 ^ mm^/sec. 



(20) 



(21) 



Therefore, the correlation-related term given in (|20p may 
be larger than (or comparable to, at least) the coefficient 
7(T under the present physical conditions. This means 
that the term should be dominant in rate equation ([1^ 
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thus totally changing the time evolution of cells from the 
case of rigid membranes. An example to take note of is a 
situation wherein Sq and 6Scui- > arc sufficiently large 
to make the term — jpoaC dominant in Eq. (|19p . In this 
case, we obtain SSdn < regardless of the values of n 
or 5*0, which may result in the cell being stationary, i.e., 
SSdii + SScur = 0. Therefore, one could stabilize the cell 
of a certain area 5*0 by imposing appropriate curvature 
increments SK whose value mainly depends on 5*0 and 7. 
This stability condition is totally different from that of 
cells confined in rigid curved membranes where neither 
5*0 nor 7 but the sign of K and the value of n are relevant. 













/r \ 







V. SUMMARY AND PERSPECTIVES 

Our results are based on the assumptions that all cell 
boundaries have invariant material constants 7 and a, 
and the shape of each jth boundary curve is described 
by a constant kj. In realistic foam, these parameters 
are determined by the nature of liquid films that retain 
three-dimensional geometry across the gap containing the 
foam. In other words, each cell boundary is not a truly 
one-dimensional curve, but a three-dimensional film with 
finite thickness whose value varies spatially within the 
film [41I, [43, H] ■ Therefore, membrane deformation 
will induce spatial fluctuations of 7 and a over the foam 
and that of Kj in each jth film. The consideration of 
these fluctuations may enhance the quantitative precision 
of the coarsening theory we have developed. 

It is interesting to point out that membrane deforma- 
tion may cause a flow of liquid through the films, as anal- 
ogous to foam drainage in response to gravity and capil- 
larity [ll|, [sO, IMI, IHj H^l ■ In fact, the deformation of the 
confining membranes leads to the rearrangement of the 
liquid film network as well as cell configuration, and thus, 
inducing pressure gradient in the liquid. Recently, it 
was shown that fluid flow on a surface with time- varying 
surface curvature exhibits three kinds of dynamical re- 
sponses depending on geometric and material constants 
[54| . In this context, we conjecture that deformation- 
induced fluid flow in the current system behaves differ- 
ently from the case of a rigid membrane, which gives rise 
to three-cornered coupling of surface deformation with 
gas transfer and associated fluid flow. Further examina- 
tion incorporated with the theory presented in Ref. (s^ l 
will shed light on the issue. 

In conclusion, wc have studied the effect of surface 
deformation on the coarsening dynamics of foam con- 
strained in a gap of flexible membranes. The growth 
rate equation of cells has been formulated by taking into 
account the correlation between surface deformation and 
gas transfer between adjacent cells. We have found that 
the correlation-related term should be dominant in the 
resulting rate equation under realistic conditions, imply- 
ing the possibility of artificial control of the coarsening 
dynamics through the confining membrane deformation. 
It is hoped that the physical picture considered in this 



FIG. 1: Illustration for calculating the Gaussian curvature K 
of a spherical surface with radius R. The north pole P is 
enclosed by a geodesic circle with radius r, i.e., by the locus 
of all points whose geodesic distances from P equal to r. The 
value of K on P is defined by Eq. (|XT|). 



article can provide a starting point for the analysis of the 
coupling between the membrane deformation and coars- 
ening dynamics in two-dimensional cellular structures. 
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APPENDIX A: DEFORMATION-INDUCED 
INCREMENT OF CELL AREA 

In this Appendix, we derive Eq. which is the ex- 
pression of the deformation-induced increment ^S'cur of 
a cell area. For a general curved surface, the Gaussian 
curvature of a point on the surface is defined by 



K = lim 



3 [27rr - C{r)] 



(Al) 



where r is the radius of a geodesic circle around the point 
and C{r) is the Icnght of its perimeter. For a spherical 
surface with radius R, for instance, it follows from Fig. [T] 
that C{r) = 27ra, a = Rs'md and r = Rd. Therefore, we 
obtain K = 1 /i?^ > and 



Kr 



(A2) 



Equation ((A2|) holds not only when if > but if < 0, 
as far as K is constant with the circular region. 
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Now we consider the area S of the circular region. It 
is given hy S = C{r' , K)dr' , and thus 



1 — cos ( VKr 



K 



Ehminating r from Eqs. (|A2p and (|A3|) yields 



SiK,C) 



2tt ~ V47r2 - K£^ 
K ' 



(A3) 



(A4) 



which converges to C'^ /{Att) in the limit of K 0. Fi- 
nally, we obtain 



dS 
dK 



2tt V47r2 - /\ £2 



£2 



X2 



and equivalently, 



2i^V47r2 - A'/:2 

(A5) 



95 



(A6) 



£ 2(27r- ifS*)' 
which completes the proof. 



APPENDIX B: SUCCESSIVE RELATION FOR 
INTERNAL PRESSURE 



Eq. dUD as 



SS, 



Sn-l 



(Bl) 



Next, Spdii is derived by substituting Eq. ([TH)) into the 
Gauss-Bonnet theorem, which leads to 



where = £j — ^q. The summations in Eq. (|B2p are 
canceled out if the sign of the summed terms is positive 
or negative depending on j. As a result, we obtain 



Spdif = -cr 



(B3) 



The internal pressure at t = ^St (/i > 1) is given which gives = +(5pdif +<5pcur for all /x > 1 through 
by the following procedure. First, Spcui- is deduced from the successive relations. 
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